- 1 

I 


AD/  A - 000  025 

AN  INTERFEROMETRIC  INVESTIGATION  OF 
SHOCK  STRUCTURE  AND  ITS  INDUCED 
SHOCK-TUBE  BOUNDARY  LAYER  IN  IONIZED 
ARGON 

Peter  I,  Brimelo,  et  al 


Toronto  University 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY 


UNCLASSIFIED 


-SECURITY  CLASSIFICATION  C,-  ‘Hli  PAGE  'IIT,,n  Putt  Fntrrr-tl 


| REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COVI’I-KTIST,  FORM 

rw-'i  74-J752 

2 GOVT  / -CEbSlON  NO- 

1.  RECIPIENT'S  catalog  l.jMdER 

fiiyrfdoc  oa5~ 

*.  TITLE  f«m/  Submit) 

AH  INTERFEROMETRIC  INVESTIGATION  OF  SHOCK 
STRUCTURE  AND  ITS  INDUCED  SHOCK- TUBE  BOUNDARY 

' 

INTERIM 

LAYER  IN  IONIZED  ARGON 

t.  PERFORMING  0*1.  REPORT  NUMBER 

T.  AUTHOR(»> 

4.  CON^RAtf  OR  OPAn'V  NUMSfcRfi) 

PITER  I BRIMELO 
1 I GLASS 

AFOSR  72-2274 

•.  performing  organization  name  ano  AODNIII 
UNIVERSITY  OF  TORONTO 
INSTITUTE  FOR  AEROSPACE  STUDIES 
ONTARIO,  CANADA 

10.  PROGRAM  ELEMENT,  PROJECT,  T ASK 
AREA  I WORK  UNIT  NUMBERS 

681307 

9783-06 

61102F 

tt.  CONTROLLING  OFFICE  NAMC  ANO  AOOREM 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NA 

It.  REPORT  DATE 

July  1974 

J.4UU  WILSON  BOULEVAKD 
| ARLINGTON , VIRGINIA  222-0 

l«  MI1MRSR  OF  PAOEt 

1 u.  MONITORING  AOKNCY  NAM  K A AQORltlirU  dttfrffit  fwn  Com  ft  Uni  OWf) 

II.  (ECuRlTV  CLAM,  f •!  Hill  nyatl) 

UNCLASSIFIED 

ik«,  auificationToownoraoino 

It.  DISTRIBUTION  STATEMENT  (•!  (Ala  ktptit) 

Approved  for  public  release;  distribution  unlimited 

• 

IT.  DISTRIBUTION  IT  ATtMCNT  fa/  (ha  aka/rant  anlarari  In  Slack  20,  ll  Mlltttnl  Ifam  BaparlJ 

It.  SUPPL  IMCNTANV  NOTH 

SHOCK- WAVE  STRUCTURE  SHOCK- BOUNDARY  LAYER  INTERACTIONS 

INTERFEROMETRY  m.,, „v 

INDUCED  BOUNDARY  LAYERS  NATIONAL  TECHNICAL 

SHOCK- TUBE  FLOW  INFORMATION  SERVICE 

IONIZED- ARGON  FLOW  " * 

20.  AatTNACT  fCanflnua  an  ravaraa  aid*  II  naeaaaaiy  and  Htntlly  by  aleak  numbarj 

The  structure  of  ionizing  shock  waves  and  their  induced  boundary  layers  in  argon 
have  been  studied  experimentally  in  the  UTIAS  4"x7"  hyporvelodty  shock  tube 
using  dual-wavelength  Interferometry.  The  relaxation  lengths  and  transitions  of 
shocks  of  Mach  numbers  16.5  and  13.5  for  initial  gas  pressures  of  about  5 torr 
have  been  measured  at  2cm  and  8cm,  respectively.  Thu  transition  profiles  and 
relaxation  lengths  are  in  reasonably  good  agreement  with  analysis.  The  effect 
of  0.47.  hydrogen  addition  to  the  argon  test  gas  for  the  stronger  shock  cusc  is 
a reduced  relaxation  length  from  20mm  to  6mm,  which  can  be  explained  on  the  basis 

DD  ,'$*“»  1473  COITION  OF  1 NOV  ft  1*  OBSOLETE  J UNCLASSIFIED 


IFCIIRIT  V CLASSIFICATION  OF  THIS  P AOE  CW*n  Cam  tn  e«rm 


UNCLASSIFIED 


liteufcrt'V  CL  AiSiriCATlCN  o*  This  I-  G£flVJi*n  C«l«  F.nlefd) 

— 


of  collision  theory.  In  the  nonequt Librium  relaxation  region  a much  faster 
ionization  rate  than  in  the  free  stream  was  found  near  the  wall.  This  phenomonor 
has  not  been  explained.  Boundary  layer  density  and  degree  of  ionization  profile^ 
at  4cm  and  20cm  behind  a shock  wave  of  Mach  number  16.5  indicate  a laminar 
boundary  layer  growing  from  0.7mm  to  2mm  in  thickness.  Similar  profiles  are 
presented  for  e station  10cm  behind  a shock  wave  of  Mach  number  13.5.  The 
profiles  are  in  better  agreement  with  analysia  for  frozen  or  nonequilibrium 
boundary  layers  than  for  equilibrium  profiles. 


itU*  UNCLASSIFIED 


AN  INTERFEROMETRIC  INVESTIGATION  OF  SHOCK 


STRUCTURE  AND  ITS  INDUCED  SHOCK-TUBE 
BOUNDARY  LAYER  IN  IONIZED  ARGON 

*y 

Pater  I.  Brimalow 


July,  197k 


ITT  IAS  Technical  Note  No.  167 
CN  ISBN  0082-5263 


Acknowledgements 


The  guidance  and  supervision  afforded  by  Dr.  I.  I.  Class  in  tho  couiv. 
of  my  Masters  degree  program  is  very  gratefully  acknowledged. 

I wish  to  extend  sincere  thanks  to  Mr*  B.  T.  Whitten  whose  helpful 
advice  and  assistance  in  every  aspect  of  this  research  project  was  Invaluable 
in  its  successful  conpletion. 

The  technical  assistance  given  by  Mr.  P.  Crouse  during  the  experimental 
program  is  greatly  appreciated. 

I should  like  to  thank  Dr.  W.  S. Liu  for  providing  the  theoretical  re- 
sults presented  in  this  report. 

The  financial  assistance  received  from  the  Air  Force  Office  of 
Scientific  Research  under  Grant  No.  M-AFOSR- 72-22?*®  and  from  the  National 
Research  Council  of  Canada  is  acknowledged  with  thanks. 


il 


Summary 


The  structure  of  ionizing  shook  waves  and  their  induced  boundary  layers 
in  argon  have  been  studied  experimentally  in  the  UIL&S  V x 7"  hypervelocity  sho-k 
tube  using  dual-wavelength  interferometry. 

The  relaxation  lengths  aid  transitions  of  ahocks  of  Mach  numbers  16.5 
and  13*5  for  initial  gas  pressures  of  about  5 to rr  have  been  measured  at  2 cm  and 
8 cm,  respectively.  The  transition  profiles  and  relaxation  lengths  are  in  reason- 
ably good  agreement  with  analysis. 

The  effect  of  0.4^  hydrogen  addition  to  the  argon  test  gas  for  the 
stronger  shook  case  is  a reduced  relaxation  length  from  20  ran  to  6 mm,  which  can 
be  explained  on  the  basis  of  collision  theory. 

In  the  nonequillbriua  relaxation  region  a much  faster  ionization  rate 
than  in  the  free-  stream  was  found  near  the  wall.  This  phenomenon  has  not  been 
explained. 


Boundary  layer  density  and  degree  of  ionization  profiles  at  b cm  and 
20  cm  bahind  a shook  wave  of  Mach  number  16.5  indicate  a laminar  boundary  layer 
growing  from  0.7  om  to  2 mm  in  thickness.  Similar  profiles  are  presented  for  a 
station  10  cm  behind  a shock  wave  of  Mach  numbor  13.5* 

The  profiles  are  in  better  agreement  with  analysis  for  frozen  or  non- 
equilibrium boundary  layers  than  for  equilibrium  profiles . 
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A 


Notation 


a speed  of  sound 

A shock  tube  cross-sectional  area,  argon  atom 

0^  specific  heat  (constant  pressure) 

Cv  specific  heat  (constant  volume) 

D A attbipolar  diffusion  coefficient 

hj  ion  enthalpy 

h enthalpy 

H total  enthalpy 

k Boltzmann  constant 

K thermal  conductivity 

Le  Lewis  number 

m mass 

M Mach  number 

n particle  number  density 

p,  P pressure 

Pr  Prandtl  number 

Q collision  cross-section 

R gas  constant 

Sc  Schmidt  number 

t time 

u particle  velocity  in  x-direction 

ua  shock  velocity 

v particle  velocity  in  y-direction 

w species  net  production  rate 

x streamwiae  direction 

y co-ordinate  direction  perpendicular  to  wall 

a degree  of  ionization 


v 


i.  ii.trod’jctio:; 


In  recent  years  much  effort  has  been  devoted  to  the  study  of  hot  g-.r 
flows,  where  the  presence  of  atoms,  ions  and  electrons  as  separate  species  in 
the  gas  considerably  complicates  the  physical  mechanisms  involved  in  the  flow. 
Therefore,  et\  understanding  of  the  behaviour  of  the  different  species  and  th^ir 
interactions  with  each  other  has  become  essential  in  predicting  the  fluid  prop- 
erties of  hot  gas  flows. 

The  relevance  of  such  research  to  present-day  requirements  is  ar.  groat 
as  ever  with  the  development  of  such  vehicles  as  the  space  shuttle.  A knowledge 
of  the  physical  characteristics  of  dissociated  and  ionized  gas  flows  is  impor- 
tant not  only  for  the  design  of  such  vehicles  that  must  withstand  re-entry  con- 
ditions, where  temperatures  caused  by  shock  and  viscous  heating  are  of  sufficient- 
magnitude  to  make  dissociation  and  ionization  effects  significant,  but  also  for 
the  test  facilities  that  are  used  to  simulate  suitable  flight  conditions  for 
these  spacecraft. 

Thu  shock  tube  istxne  such  facility  which  is  capable  of  generating  high 
temperature,  high  Mach  number  gas  flows  in  which  dissociation  and  ionization 
processes  take  place.  The  UTIAS  U"  x 7"  hypervelocity  shock  tube,  which  vac 
used  for  this  experimental  study,  is  capable  of  very  high  performance,  generating 
shock  Mach  numbers  in  excess  of  24  in  argon  gas.  The  relatively  large  cross- 
sectional  area  of  this  shock  tube  makes  it  especially  suitable  for  the  study 
of  boundary  layers,  the  effect  of  which  on  the  free-stresm  is  virtually  negligi'l 

Argon  is  chosen  as  a test  gas  for  a number  of  reasons.  It  is  a mon- 
atomic gas  and  consequently  dissociation  effects  are  eliminated,  which  simpLLfi 
the  experimental  technique,  as  will  be  seen  later.  Argon  has  a relatively  low 
first  ionization  potential  (15.7  eV)  compared  to  some  other  monatomic  gases,  whi 
practically  speaking  means  that  for  a given  set  of  initial  conditions  a higher 
degree  of  ionization  is  achieved  than  for  a gas  of  higher  potential.  The  thermo- 
dynamic and  transport  properties  of  argon  have  been  well  established  over  a via: 
range  of  tenperaturos  and  pressures.  Lastly,  argon  1q  available  commercially  and 
relatively  cheaply  with  a high  degree  of  purity  and  it  is  a safe  gas  to  work  with. 

The  experimental  study  described  in  this  report  was  originally  intend > 
to  investigate  only  the  shock- tube  wall  boundary  layer  formed  behind  ionizing 
shock  waves.  However,  it  soon  became  apparent  that  the  optical  system  which  vu 
set  up  to  study  the  thin  boundary  layer  region  was  equally  suited  to  studying 
the  shock-wave  structure  with  which  it  initially  interacts. 

In  the  literature  many  theoretical  analyses  have  been  presented  describ- 
ing the  structure  of  ionizing  shock  waves  (e.g.,  Hefs.  1-4)  and  also  some  on 
none tationary  ionized  boundary  layers  (Refs.  5-7)  but  very  little  has  beer,  re- 
ported on  experimental  results  obtained  in  the  areas  of  interest  here,  Rel'erer. 1 
4,  8 and  9 give  experimental  results  obtained  for  the  boundary  layer  region 
while  Refs.  3 and  10  give  some  experimental  profiles  for  the  shock  wave  radar.:..1  ! - 
region. 


1 


ANALYSIS 


The  theoretical  analyses  presented  in  this  section  are  not  meant  to  be 

complete  treatments  of  the  problems,  but  more  as  an  outline  of  the  relevant  equa- 
tions and  factors  which  are  important  in  obtaining  a general  solution. 

The  flow  regions  may  be  separated  into  distinct  sections  as  follows 
(Fig.  l).  In  the  freestream  the  Rankine-Hugoniot  diock  front,  behind  which  the 
flow  is  frozen,  is  followed  by  none q'd librium  flow  during  which  energy  equilibria- 
tion  (or  relaxation)  processes  take  place.  The  chemical  kinetics  involved  in  this 
region  are  discussed  later.  Behind  the  relaxation  region  the  gaB  is  in  thermo- 
chemical  equilibrium,  although  some  changes  do  take  place  slowly  due  to  radiation 
losses.  The  viscous  boundary  layer  growth  begins  at  the  shock  front  and  extends 
back  through  the  nonequilibrium  flow  region  into  the  equilibrium  flow  region.  The 
theoretical  boundary  layer  analysis  presented  in  Sec,  2,2  is  limited  to  the 
equilibrium  boundary  layer  problem. 

2.1  Freestream  Relaxation  Region 

The  governing  equations  of  mass,  momentum  and  energy  written  for  one- 
dimensional  steady  state  flow  in  the  shock- wave  region  are  (Refs,  11  and  12) 

•ft  (PU)  - 0 (2.1.1) 

(pu2  + p - 3 u ) - 0 (2.1.2) 

4C  . _ du  4 / du  Y+  L.  (y  £ 'N\  to  i 

P dx  P dx  3 ^ \ dx  / dx  \ dx  / (2.1.3) 

The  momentum  equation  (Eq.  2.1,2)  is  derived  from  the  Navi er- Stokes 
aquations.  The  4/3  u coefficient  arises  from  Stokes'  assumption  that  the  second 
coefficient  of  viscosity  \ « -2n/3  and  noting  that  for  a monatomic  gas  the  co- 
efficient of  bulk  viscosity  is  zero, 

Tne  terms  on  the  right-hand  side  of  the  energy  equation  (Eq.  2,1.3) 
describe  entropy  production  in  the  shock  wave. 


Equations  2.1.1,  2,1.2  and  2.1.3  describe  the  flow  both  through  the  shock 
front  (translational  excitement  only)  and  through  the  ensuing  relaxation  region. 

In  the  former  region  (layer)  the  gradients  in  tho  flow  variables  are  large  and 
the  effects  of  viscosity  and  heat  conduction  must  be  considered.  In  this  analysis 
we  shall  be  concerned  only  with  the  relaxation  region  and  since  the  gradients  in 
flow  properties  are  not  large,  the  viscous  and  heat  conduction  dissipative 
processes  may  be  neglected.  Thus  Eqs.  2.1.1,  2.1.2  and  2.1.3  may  be  reduced  to 
taka  the  following  form* 


d_ 

dx 


(P^) 
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(2.1.4) 


pu 


d£ 

dx 


(2.1.5) 
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0 


(2.1.6) 


dx  dx 

which  when  integrated  yield  the  standard  equations  (Ref.  13) 


pu  - PjUj^ 

(2.1.7) 

p + pu2  - + P1U2 

(2.1.8) 

\ 

h + i u2  ■ h^  + i u2 

(2.1.9) 

where  u and  axe  measured  relative  to  the  fixed  shock  front. 

Following  Ref.  1,  Eqs.  2,1.?,  2,1,8  and  2,1,9  together  with  the  two  - 
temperature  concept  in  the  equation  of  state  (for  the  electrons  and  heavy 
particles) 


* 


p - pR  (T  + a T ) (2.1.10) 

and  enthalpy  equation 

h - 5/2  R (T  + a l'e)  + a R 8j  (2.1,11) 

are  solved  to  give  the  particle  velocity  (assumed  the  same  for  heavy  particle  and 
electron  alike) 


(2.1.12) 


where  the  negative  sign  gives  the  physically  meaningful  result  and  Mg,  the  shock 
Mach  number  is  defined  as  * 


ua  / (5/3  R T^1/2  - ug/a1 


The  gas  tenperature  T is 

T 


M2 


-HvCO  J) 


(2.1.13) 


To  determine  the  degree  of  ionization,  a and  electron  temperature  T , 
we  now  consider  the  mass  and  energy  conservation  equations  for  the  electron  gas 
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L 


- 1 


r 


5-  (n  u)  = ft 

dx  ' e e 


nu  i 13/2  k T ) + n k T 
e dx  <1  e e dx 


(2. 1.1*0 

3-e  ( r)  *VeA+  * 3/2  V?ri 


where  x is  the  atreamwise  co-ordinate  and  v .,v  _ are  the  elastic  collision  fre- 

eA  el 

quencies  for  electron-atom  and  electron-ion  collisions  respectively  and  are  related 
to  the  corresponding  elastic  collision  cross-sections  as  (Ref.  2) 


eA 


**(i-q)p  ( 8 k Te  \ 
3 \ vm  J 


1/2 


Q, 


, T • 22 

el  m„ 


(BkTe-) 

V / 


e 

1/2 


’eA 


(2.1.15) 


*el 


Equations  (2,1,14)  may  be  written  in  terms  of  the  degree  of  ionization  a,  as 


d a i ' * 

u s \ * a. 


dT  a,, 

“ 3T  * *.  S 


3(^)(W<*-T.>-  ET6] 


(2.1.16) 


where  a ^ and  a0  are  the  ionization  rates  due  to  atom-atom  and  atom- electron  colli- 
sions, respectively.  Then 


da  • , 

e 


and 


dl  /m\V  +v_ 

2 - 2 ( -S  ) -2A Si  (t-T  ) 

dx  V mA  J u u V 


2 ^ fi 

3 a u 


(2.1.17) 


(2.1.18) 


If  the  local  steady  state  approximation  is  assumed,  i.e.,  dT  /dx  » 0, 
then  Eq.  2.1,18  gives  " 


, .T  I ( V \ <V  3/2  V 

e 3 \ / »oA  * v8l 


a \ 
cT  ) 


(2.1.19) 


A 


wherci  a.  , v . and  v _ are  all  functions  of  T (see  Ref.  1 for  a ). 
e’  eA  el  e ' e 

In  principle,  Eqs.  2.1.17  and  2,1.18  or  Eqs.  2.1.17  and  2.1.19  may  be 
olved  numerically  to  give  a and  T for  given  x,  although  the  boundary  condition!’ 

for  Tg  require  careful  consideratiSn. 

The  chemical  kinetic  theory  describing  the  processes  in  a shock  wave 
leading  to  excitation  and  ionization  has  been  discussed  by  a number  of  authors 
(e.g.,  Refs.  1,3,4).  A sinple  outline  is  given  here  of  the  more  popular  theory 
which  involves  a two-step  atom- at am  collision  process  followed  by  a two-step 
clcctroi-atom  and  single-step  electron- electron  collisional  model  (see  Fig.  l). 

With  the  initial  discontinuous  jump  in  flow  conditions  having  taken 
place  through  the  Rankine-Hugoniot  shock  front  (whose  thickness  is  of  the  order 
of  a few  mean  free  path  lengths),  the  gas  is  in  a frozen  state,  the  energy  trans- 
fer taking  place  only  in  the  translational  mode.  The  relaxation  region  begins 
with  a two-srep  atom-atom  collisional  process 

A + A -+  A + A* 

A*  + A -»  A + A++e 

wJdch  takes  some  atoms  to  their  first  excited  state  before  further  collisions 
energize  them  sufficiently  to  free  the  electron.  This  relatively  slow,  in- 
efficient process  leads  to  the  production  of  a small  number  of  electrons. 

The  next  phase  is  characterized  by  a two-step  and  a single-step  event 
by  which  the  energy  transfer  rate  is  considerably  speeded  up. 

A + e # A*  + e 
A*  + e A+  + e + e 

A + e A + e + e 


The  latter  process  soon  predominates,  leading  to  the  rapid  producticn 
of  electrons  until  an  equilibrium  state  is  reached  with  electrons  recombining 
at  the  same  rate  at  which  they  are  produced.  This  balance  is  somewhat  disturbed 
by  radiative  cooling  by  which  some  energy  is  lost  through  the  amission  of  radiation 
with  a consequent  drop  in  electron  number  density.  This  process  affects  the  entire 
flow  uniformity  behind  the  shock  front. 

2.2  Equilibrium  Boundary  Layer 

The  boundary  layer  formed  behind  a normal  shock  wave  travelling  along 
an  infinite  flat  plate,  as  is  the  case  for  the  shock  tube  boundary  layer,  is 
nonstationary. 

The  nonstationary  laminar  boundary  layer  equations  for  a partially  ionize 
gas  have  been  treated  by  many  authors  and  may  be  written  as  follows  (Ref.  6) 

If  + lx  (pU)  + |y  (pv)  “ 0 (2*2.1) 
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; } 


i 
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whe.  • 


and 


( bu  , bu  , bu  N 

p K’St  + u v b£  ) 

t? ) 

(2.2.2) 

- 0 

(2.2.3) 

nfb H . „ bn  . „ bH  N 
PVbt  + u b£  + v by. 

" b [ ’ Vqd  + 

“ &(r). 

(2.2.4) 

q_c  ■ -K  ^ , the  conductive  heat  flux 

qd  " “Pl>A  (Qj.  + 5/2  RT)  ^ , the  diffusive  energy  flux,  assuming 
amblpolar  diffusion. 


Equations  (2.2, 1-2. 2.4)  are  respectively  the  continuity,  x and  y momentum 
and  energy  equations. 


These  equations  may  he  reduced  to  steady  state  equations  (b/bt  «0)  by  a 
transformation  of  the  co-ordinate  system  and  changing  the  boundary  conditions,  as 
shown  in  Fig.  2.  Equation  2,2.3  is  a statement  of  the  fact  that  the  pressure  is 
invariant  through  the  boundary  layer,  (terms  of  the  second  order  being  neglected). 

If  we  assume  that  the  freestream  pressure  is  independent  of  x,  then  the  steady- state- 
equations  are 


a 

b£ 

(PU)  + (PV)  -0 

(2.2.5) 

bu 

bu  b f bu  \ 

(2.2.6) 

pu  s* 

PVby“  by  by,) 

- 1 **  1 

b [ b Ai2  \ 
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If  the  gas  remains  in  theraochemical  equilibrium,  the  Saha  equation 
(Eq.Al)  can  be  used  tb  directly  relate  the  degree  of  ionization  a to  tenperature 
T,  for  a given  pressure.  In  this  case  the  energy  equation  (Eq.  2,2 »7)  may  be  re- 
duced to  the  familiar  classical  form  (Ref,  l4) 
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provided  that  the  Prandtl  number,  Pr  relating  viscosity  and  heat  conduction,  is 
defined  as 


Pr 


Models  for  the  transport  properties 
suits  obtained  there  we  have 


are  developed  in  Appendix  A* 


Using  the  re- 
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(2.2.10) 
(2.2.11) 


Note  that  the  definition  of  the  total  thermal  conductivity  includes  the  normal 
thermal  conductivity  K,  as  well  as  a diffusive  energy  transport  term. 

In  classical  boundary  layer  theory, Pr  is  generally  assigned  a constant 
value  since  it  is  virtually  Independent  of  temperature  at  the  usual  temperatures 
considered.  However,  in  the  case  of  ionised  boundary  layers  the  range  of  teiqpera- 
ture from  the  cold  wall  to  the  hot  frees tream  is  so  great  that  variations  in  vis- 
cosity, specific  heat  and  thermal  conductivity  are  such  as  to  cause  significant 
changes  in  Pr. 

If  the  assuaption  of  thermoohealcal  equilibrium  is  not  made,  another 
formulation  of  the  energy  equation  may  be  obtained  (Ref.  6)  i.e., 

pu  S*PT  3y  " 35  (fe'  3y)  + 35  [ “ O' ' f ) 35(f)] 

* 35  [ fc'  (I-1)  hi  35  ] (2-s-u) 
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whore  Pr'  is  based  on  values  of  C'  and  K'  defined  as 

P 


<r  - 5/2  r (1  + a) 

K'  - K 

In  Eq.  2.2.12,  the  Levis  number  Is  defined  as 

P DAPr' 


(2.2.13) 

(2.2.14) 


(2.2.15) 


where  I>A  is  the  asibipolar  diffusion  coefficient.  The  ratio  of  Prandtl  number  to 
Levis  number  is  referred  to  as  the  Schmidt  number,  Sc 


i.e., 


Sc  “ fr  “ pfc 


(2.2.16) 


In  addition  to  the  extra  term  which  appears  in  the  energy  equation,  an 
equation  for  species  conservation  must  be  solved.  For  ions,  this  is 


pu  i!  + pv  3£*ly(pDA  ) + * (2.2.17) 

assuming  ambipolar  diffusion,  where  w is  the  speaies  net  production  rate.  For 
equilibrium  flow,  of  course,  Eq.  2.2.17  may  be  replaced  by  the  Saha  equilibrium 
coJQOsltlon  relation  (Eq.  A.l).  However,  for  nonequilibrium  solutions,  Eq.  2.2.17 
must  be  considered  along  with  the  energy  equation  2.2.12. 

The  theoretical  results  presented  for  the  frees tream  relaxation  region 
were  obtained  by  Liu,  based  on  an  extension  of  the  analysis  presented  in  Sec.  2.1. 

A full  description  of  this  is  presented  in  Ref.  15. 

The  theoretical  equilibrium  boundary  layer  results  given  in  this  report 
were  obtained  by  Whitten  (Ref.  16),  while  a coiqplete  boundary  layer  analysis  for 
equilibrium,  frozen  and  nonequilibrium  flow  is  presented  in  Ref.  15. 

3.  DESCRIPTION  OF  EXPERIMENTAL  EQUIPMENT  AND  TECHNIQUES 

Since  the  original  construction  of  the  UTIAS  4"  x 7"  hypervelocity  shock 
tuba,  several  improvements  have  been  made  to  the  facility,  increasing  both  its 
performance  and  capability  to  the  level  required  for  the  successful  outcome  of  the 
experimental  program  described  in  this  report.  These  modifications  concern  the 
driven  section  vacuum  system  and  the  electronic  timing  circuits  associated  with  shock 
speed  measurement  and  triggering  of  the  optical  system's  light  source.  Figure  3 
illustrates  the  test  section  end  of  the  facility. 

3.1  The  Shock  Tube 

The  design  and  construction  of  the  UTIAS  4"  x 7"  chock  tube  is  described 


\ 
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in  Ref.  17. 

A combustion  driver  using  a stoichiometric  mixture  of  7.%  (by  pressure) 
oxygen,  20&  hydrogen  and  72,%  helium,  was  used  to  create,  by  combustion,  the  re- 
quired driver  pressure  P^.  This  mixture  was  determined  by  Bristow  (Ref.  18)  to 
give  a smooth  combustion4 and  an  acceptable  rise  time  and  peak  pressure.  The  gas 
mixture  was  ignited  impulsively  by  means  of  a tungsten  wire  (0.015"  dia.)  located 
down  the  centre  of  the  driver  and  across  which  a 13  Kv  potential  was  applied. 

The  ratio  of  combustion  pressure  to  test  gas  pressure  P^  was  found  by 
Bristow  to  be  related  to  the  incident  ahook  wave  Mach  number  by  the  empirical 
expression,  Ln(PjVP.)  ■ [0.1*22  Mfl+  4.02]  in  the  range  10  < Mg  < 24.  The  factor 
by  which  the  stoichiometric  gas  mixture  pressure,  P._  is  increased  by  combustion 
is  approximately  6.7. 

Design  of  the  dlaphragaa  was  an  important  factor  in  the  series  of  experi- 
ments conducted  here  since  repetition  of  flow  conditions  was  required  in  order 
that  the  boundary  layer  could  be  studied  at  different  stations  behind  a shock  wave 
of  given  strength.  Accurate  prediction  of  the  shock  speed  was  also  important, 
as  will  be  seen  later,  in  triggering  the  light  source  for  the  optical  equipment. 

The  diaphragms,  made  from  type  304  stainless  steel  and  annealed,  were 
scribed  to  ensure  proper  petaling  when  burst.  Poliowing  Brlatow,  whose  calibra- 
tion of  this  type  of  diaphragm  ia  given  in  Ref.  l8,  the  bursting  pressure  P^  is 
related  to  the  diaphragm  thickness  t,  the  thickneaa  after  scribing  h,  the  un- 
supported diameter  d and  the  ultimate  tensile  strength  of  the  material  a,  by 
the  equation 


p4 


where  K is  a factor  close  to  unity,  depending  on  t and  n ia  an  experimentally  de- 
termined value  (n  2.2).  (see  Ref.  l8  for  specific  values  of  K). 

The  driven  section  of  the  tube  is  48  ft  in  length,  the  test  section 
being  located  approximately  4$  ft  from  the  diaphragm  station.  The  material  used 
for  construction  of  the  driven  section  was  a high  strength  alloy  steel  (type 
SPS-245)  with  the  Internal  surfaces  being  ground  to  a 32  micro-inch  finish  and 
plated  with  a 0.002  inch  layer  of  chromium. 

Between  experiments  the  tube  was  blown  out  with  high  pressure  air  and 
cleaned  with  acetone  and  trichlorethylene  so  as  to  maintain  a reasonable  degree 
of  cleanliness. 

The  vacuum  pump  system  for  the  driven  section  of  the  tube  consists  of  a 
Kinney  type  K3-47  mechanical  pump,  a Hereaus  Rootea  type  VP-RQ-350A  pump  and  a 
CVC  type  PMCS-6B  diffusion  pump  which  in  series  were  capable  of  bringing  the  tube  pressure 
down  to  10-°  nmHg,  with  an  outgassing  rate  of  about  1*10“*  nmHg/min,  Generally,  the 
tube  was  pumped  down  for  3 to  4 days  between  experiments  to  ensure  as  good  a vacuum 
as  possible. 

The  argon  test  gas  (with  a purity  of  99,99%  by  volume)  was  introduced 
directly  into  the  driven  section  immediately  before  an  experiment waa  carried  out, 
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and  the  pressure  monitored  by  a Wallace -Tier nan  gauge  (0-50  mnHg  range).  In  each 
experiment  a test  gas  initial  pressure  of  5 moHg  was  set  by  this  gauge.  A more 
accurate  measurement  of  the  pressure  waB  taken  using  an  oil  manometer,  isolated 
before  the  experiment  was  conducted.  This  instrument  is  described  in  Ref.  18, 
along  with  its  ' calibration. 

For  the  one  experiment  in  which  it  was  required  to  add  a small  per- 
centage (0.4#  by  pressure)  of  hydrogen  to  the  test  gas,  a reservoir  of  known 
volume  was  required.  The  following  procedure  was  adopted.  Using  the  known  volumes 
of  the  dump  tank  (35.4  ft3)  a volume  Vg  is  calculated  as  1*5 .0  ft3.  With  an  arbi- 
trarily set  pressure  of  40  rnaHg  (P^)  in  the  hydrogen  gas  line  (as  measured  by  a 
Wallace  & Tieman  gauge),  the  hydrogen  was  released  into  the  tube  and  the  pressure 
(Pg)  measured  with  a Mcleod  gauge. 

Using  the  standard  equation,  PjV^  ■ PgVg,  the  volume  (V^)  of  the  hydrogen 
gas  line  (including  pressure  gauge  volume)  was  calculated. 

3.2  Shock  Velocity  Measurement 

Measurements  of  shock  velocity  were  obtained  between  a reference  station 
and  4 positions  in  the  tube  using  standard  piezoelectric  pressure  transducers. whose 
outputs  were  amplified  before  being  applied  to  2 Hewlett-Packard  (type  3734  A)  and 
2 Rao&l  electronic  counters  (see  Fig.  4),  measuring  in  micro-second  units*  A fifth 
velocity  measurement  was  teken  across  two  feet  of  the  test  section  using  a Hewlett- 
Packard  (type  5325  A)  counter  measuring  in  0.1  micro-second  units. 

As  mentioned  in  Ref.  18,  the  degree  of  shock  attenuation  recorded  with 
this  system  is  so  low  as  to  lie  within  the  experimental  error  of  the  velocity 
measurement  technique, 

A Kistler  pressure  transducer  was  used  (station  0,  Fig.  4)  to  record  an 
oscillogram  of  the  pressure  history  through  the  ehodk  wave  and  subsequent  gas  flow. 

3.3  Optics 

The  optical  system  employed  in  this  experimental  program  comprised  a 
pulsed  dual- frequency  laser,  a Mach-Zehnder  interferometer  and  a plate  camera. 
Figure  5 shows  a schematic  layout  of  the  optics  and  Fig.  6 gives  a plan  view  of 
the  experimental  apparatus. 

3.3.1  Light  Source 

A TRO  104a  ruby  laser,  pulsed  by  a Pockels  cell  ^-switch  and  fitted  with 
a second  hammic  generator,  provided  the  6943  A and  3471.5  A component  wavelength 
light  source.  A full  description  of  this  laser  and  its  advantages  and  disadvant- 
ages over  more  conventional  light  sources  is  given  in  Ref.  13.  The  ruby  rod  of 
this  type  of  laser  is  energized  by  a flash- Imp  for  a period  of  about  950  usees 
before  the  optimum  condition  is  obtained,  prior  to  Cl-switchlng.  A signal  from 
one  of  the  aforementioned  pressure  transducers  1b  used,  via  a delay  unit,  to 
trigger  the  flash  lamp  and  a similar  signal  from  another  tranaducsr  is  used,  again 
via  a delay  circuit,  to  trigger  the  Pockels  cell  at  a predetermined  time  (see 
Fig.  4).  In  this  manner  the  desired  region  of  flow  is  exposed  to  the  film  plates. 
This  system  of  shuttering  releases  a short  burst  of  light  in  the  order  of  20  nano- 
seconds, which  is  instantaneous  relative  to  the  shock  speed  (say  5 Km/sec).  Of 
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course  accurate  knowledge  of  tna  shock  speed  Is  required  prior  to  a run  for 
successful  operation  of  the  laser-shutter  system.  A measurement  of  the  time 
elapsed  between  the  shock-wave  arrival  at  a reference  station  and  operation  of 
the  ft- switch  is  recorded  on  a Hewlett-Packard  (type  3731*  A)  counter  (see  Fig.  4) 
and  this,  in  conjunction  with  a shock  speed  measurement  from  the  same  reference, 
enables  an  accurate  determination  of  the  flow  field  exposed  to  the  plates,  with 
respect  to  the  shock  wave*  This  Is  essential  for  cases  where  the  shock  wave  is 
outside  the  field  of  view. 

A short  focal  length  (150  an)  convex  lens,  whose  selection  is  determined 
by  the  maximum  focal  length  that  will  achieve  the  necessary  cOllimation  of  light 
to  the  required  field,  is  used  to  focus  the  laser  beam  onto  a light  diffuser  placed 
at  the  focal  point  of  the  first  parabolic  mirror  of  the  interferometer*  The 
Introduction  of  a diffuser  reduces  the  spatial  coherenoe  of  the  laser  light, 
resulting  in  better  quality  intarferograms . 

3*3.2  Mach-Zehnder  Interferometer 

The  design  and  construction  of  the  Mach-Zehnder  interferometer  used  in 
theee  experiments  is  described  in  Ref*  19.  The  adjustment  procedure  for  obtaining 
the  required  fringe  system  is  described  in  that  report.  The  alignment  procedures 
for  this  study,  involving  the  interferometer  and  the  light  source,  is  described  in 
detail  in  Appendix  B* 

3*3*3  Interferogram  Photography 

The  camera  used  for  the  photography  is  shewn  schematically  in  Fig.  5. 

It  consists  of  a shutter  ( nontrigger  able) , a beam  splitter  (60Mo£)»  a plane 
mirror  and  2 line  filters  to  separate  out  the  conponent  wavelengths.  It  is  de- 
signed to  hold  4"  x 5"  film  plates. 

Two  lenses  were  used  to  obtain  the  required  overall  magnification 
(«  1.8).  A schematic  ray  diagram  of  this  arrangement  is  shown  In  Fig.  7. 

Apart  from  the  magnification  factor,  the  lens  selection  was  made  with  several 
other  criteria  In  mind.  Thasa  may  be  aumarlaed  as  follows) 

i)  The  combination  of  the  focal  lengths  of  the  lenses  chosen  to  give 
the  desired  magnification  must  be  compatible  with  the  physical  limitations  of 
the  camera  and  space  between  the  camera  and  interferometer. 

il)  The  aperture  of  the  light  beam  at  any  station  must  fall  within 
the  bounds  of  the  lens,  filter,  splitter  and  mirror  apertures.  Working  at  the 
outer  edges  of  lenses  can  give  rise  to  fringe  distortions  and  focussing  problems. 

ill)  The  bams  must  not  be  focussed  in  the  plane  of  any  of  the  optical 
elements,  since  burning  of  the  element  may  occur. 

lv)  The  use  of  an  ultra-violet  light  source  (as  is  the  3471.5  8) 
poses  some  problems  since  this  wavelength  falls  at  the  lower  end  of  the  bandwidth 
for  which  git*  a lenses  are  useful.  The  transmission  of  light  is  poor  and  the 
thicker  the  lens,  the  more  the  light  lost, 

v)  The  overall  magnification  also  determines  the  intensity  of  light 
Available  to  the  film.  The  use  of  a second  harmonic  generator  cuts  down  on  the 
available  light  energy  from  the  laser  and  while  there  is  no  problem  with  the 
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primary  (69^3  &)  wavelength,  the  secondary  (34?1.5  S)  wavelength  requires  that 
the  laser  be  very  finely  tuned. 

Figure  8 illustrates  the  advantage  of  using  as  large  e.  magnification 
as  possible  for  the  film  plate,  rather  than  relying  completely  on  enlarging 
when  printing. 

The  film  usod  for  the  experiments  was  Kodak  Royal  X Fan  (1250  ASA) 
in  4"  x 5"  plates.  Different  developers  were  tried  (DK-50,  D-ll,  D-8)  to 
improve  the  fringe  definition  on  the  negatives.  It  was  found  that  a undiluted 
DK-50  solution  gave  the  beat  results  with  a developing  time  20$  longer  than 
that  recommended. 

Enlargements  of  between  6 and  8 times  the  negative  else  were  printed, 
giving  en  overall  magnification  of  12  to  16  over  the  laboratory  scale.  Kodak 
type  1594  grade  4 Industrial  paper  (in  rolls  of  40n  x 100')  was  used  for  this 
purpose,  Limits  on  the  size  of  prints  made,  were  lapoeed  mainly  by  the  practical 
aspects  of  developing  and  general  ease  of  handling  in  the  darkroom. 

3*3*4  Test  Section  Plane  of  Focus 

The  problem  of  fixing  a plena  of  focus  in  the  test  section  and  also 
that  of  providing  a suitably  accurate  scale  on  the  interferogrms  was  overcome 
in  the  following  manner. 

Several  holes  were  drilled  in  a metal  plate  at  approximately  l om 
intervals.  Cross-hairs  were  placed  ever  each  of  these  holes  and  the  distances 
between  cross-hairs  measured  accurately  with  a Hiiger  T.500  Universal  Measuring 
Projector  (accurate  to  0.00001"  ).  This  plate  was  fixed  in  the  interferometer 
between  the  second  mirror  and  second  splitter  (aoe  Fig.  9),  making  sure  that 
the  plans  of  the  cross-hairs  was  parpen dloular  to  the  light  beam.  The  optimum 
plane  of  focus  in  the  test  section  is  a matter  of  conjecture.  Following  Ref. 20, 
where  it  is  shown  that  refraction  errors  are  minimised  by  focussing  on  a plane 
2/3  across  the  test  section  from  the  window  nearest  to  ths  light  source,  the 
croso-hairo  wars  first  fixed  in  the  test  station  and  an  image  focus  plans  position 
established  (using  only  the  parabolic  mirror,  ?2)  before  being  moved  into  the 
ccopenaating  chamber  light  beam  path  and  adjusted  In  position  until  they  were 
again  brought  into  focus. 

The  cross-hairs  served  another  useful  purpose  in  that  of  determining 
the  true  position  of  the  shook  tube  wall.  As  can  be  seen  from  the  flow  lnter- 
ferograms,  (e.g.,  Figs.  11,12  and  13)  the  real  wall  surface  is  ill-defined  owing  to 
the  effects  of  refraction  of  the  light  through  the  boundary  layer  density  gradient. 
Before  each  experiment  the  camera  was  focussed  on  the  wall  edge  nearest  the 
camera  (which  gave  the  clearest  wall  imags)  and  ths  frlngta  blanked  out  near 
the  wall  by  means  of  a plate  inserted  in  the  compensating  arm  of  the  interfero- 
meter. A photograph  was  then  taken  (one  for  each  wavelength)  and  used  to  measure 
the  distances  from  the  cross-hairs  to  ths  wall  surface  (Fig.  10).  These  distances 
(when  adjusted  by  the  correct  magnification  faotor)  were  usod  to  determine  a wall 
position  on  the  interferogrems  (with  an  estimated  accuracy  of  + 0.05  mm),  which 
as  can  be  seen  in  Appendix  C,  i«  inportant  in  the  analysis  of  Tihe  interferogrems . 

This  method  of  using  corresponding  planes  in  the  test  section  and  com- 
pensating arms  of  ths  interferometer  is  felt  to  be  very  useful  for  experiments  in 
which  it  is  desired  that  no  obstruction  be  placed  in  the  shock  tube  that  may 
disturb  the  flow. 
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4.  RESULTS  AMD  DISCUSSION 


4.1  Shock  Front 


Measurements  of  fringe  shift  through  the  Rankine-Hugoniot  shock  front 
Vera  taken  from  the  interferograms  of  three  experiments,  i)  M ■ 16.53  (pure 
argon),  ii)  M *»  16.6Q  (argon  + 0*4$  Hg)  and  iii)  M ■ 13.59  (Sure  argon).  The 
interfe rogrsml  are  shown  in  Figs.  11,  12  and  13  respectively.  The  Rankine-Hugoniot 
front  is  visible,  on  the  right-hand  side  of  the  interf erograas , as  an  almost  dis- 
continuous shift  in  the  fringes.  The  shift  through  this  front  for  \ (about  1.5 
for  the  conditions  described  here)  is  theoretically  twice  that  for  assuming 
frozen  flow  conditions  immediately  behind  the  front  (i.e.,  a ■ 0).  Due  to  the 
very  high  density  gradient  through  the  front,  it  appears  considerably  thicker  than 
is  really  the  case.  This  vdl..  be  discussed  in  more  detail  later.  The  measurements 
of  fringe  shift  (i,e«,  th».  displacement  ot  a fringe  between  the  undisturbed  gas 
region  and  a reasonably  uniform  flow  region  behind  the  front)  were  substituted  into 
the  interfe rometrio  equations  (Eqs.  CIO  and  Cll)  and  p and  a calculated  as  shown 
in  Table  I. 


TABLE  I*  EXFER3MBBTAL  DATA 


■ 

P^CmmHg) 

S\ 

sx2 

<*($} 

5.12 

0.724 

1.520 

4.3l4»10"2 

-6 .74*10" 2 

4.38o«lo“® 

0.0 

Slsl 

5.17 

0.715 

1.599 

4.558s10‘2 

3.52-10"1 

4.4l2«10'2 

0.0 

n 

5.09 

0.744 

1.608 

4.525s10“2 

2.03S10’1 

4, 328*10" 2 

0,0 

Agreement  with  the  theoretical  density  values  p^,  as  predicted  by  the  Rankine- 
Hugoniot  relations,  is  reasonably  good  although  more  so  for  the  atronger  shook 
cases.  Verification  of  the  frozen  nature  of  the  gas  behind  the  shock  front  is 
given  by  the  negligible  degree  of  ionization* 

An  estimate  of  the  error  involved  with  these  p and  a calculations  may 
be  determined  using  the  classical  error  equation  derived  in  Ref.  21. 


“ ■ { tq}  i l * ii*{in}  'j*  ••• 


where  q Is  the  error  in  the  calculated  quantity  A,  which  is  a function  of  jhe 
measured  quantities  a^,  a2,  a^  etc. [i.e.,  A • fCa^a^a^...)]  and  etc*» 

are  the  errors  in  a^a^a^  etc. 

This  equation  may  be  applied  to  the  interferometric  equations,  Eqs.  (CIO) 
and  (Cll),  as  described  in  Ref.  22.  For  nondimens ional  fringe  shift  measurement 
errors  of  0.03,  the  errors  predicted  by  Ea.  4.1,1  in  density  p and  degree  of  ioni- 
zation a are  respectively  9»53*10"<4  Kg/nr  for  2.1$  of  the  predicted  frozen  value, 
pj  and  0.21$, 

Cm 

As  mentioned  earlier,  the  apparent  shock  front  thickness  (approximately 
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1.5  mm)  is  considerably  greater  than  the  theoretical  value.  From  Ref.  23  the  limi- 
ting value  for  shock  front  thickness  (Ms-*«)  is  1.1*2  iQ,  where  tQ  is  the  mean  free 
path  length  of  the  undisturbed  gas.  For  the  conditions  P.  «•  5 torr  and  T."  300°K, 
i = 2 x 10'14  inches  for  argon,  which  gives  a limiting  shock  front  thickness  of 
? .84  x 10"  ^ inches.  Thus  there  is  a great  difference  between  the  observed  shock 
front  thickness  and  the  theoretically  predicted  value.  This  may  be  explained  by 
a combination  of  errors.  Firstly,  considering. refraction,  the  density  gradient 
through  the  front  is  large  (typically  4,5  Kg/nr/m  for  a linear  gradient).  However 
if  we  assume  that  the  interferometer  is  perfectly  aligned  (normal  to  the  test 
section)  any  light  ray  passing  through  the  shock  front  will  be  refracted  very  quickly 
out  of  the  thin  front  region  into  the  uniform  region  behind  it.  Consequently  errors 
due  to  refraction  may  be  ignored  for  the  perfectly  aligned  interferometer.  However, 
if  the  interferometer  is  misaligned,  even  to  a small  degree,  then  refraction  effects 
become  quite  significant,  aince  many  more  rays  pass  through  the  front. 

The  alignment  procedure  described  in  Appendix  B is  not  as  accurate  about 
the  y-axls  as  It  is  about  the  x-axia  of  the  shook  tube  because  there  is  no  reference 
surface  available  on  which  to  observe  the  spurious  fringes  associated  with  misalign- 
ment, Consequently,  a two-dimensional  effect  is  inevitable  when  looking  at  the 
thin  shock  front  region. 

\ 

The  fact  that  the  shock  front  appears  considerably  thicker  for  the  weaker 
shock  may  be  the  result  of  poorer  alignmsnt  and  consequently  a larger  two-dimensional 
effect.  The  Interferometer  and  the  light  source  were  left  undisturbed  as  far  as 
possible.  However,  a major  adjustment  of  the  interferometer  was  made  prior  to  the 
weaker  shook  experiment  and  it  ia  reasonable  to  suppose  that  a poorer  alignment  was 
achieved  about  the  y-axis.  This  would  also  account  for  the  somewhat  larger  dis- 
crepancy between  calculated  and  predicted  frosen  density  values  behind  the  weaker 
shock  front. 

4.2  Frees tr earn  Relaxation  Region 

The  interferograms  shown  in  Fig.  11  illuatrata  well  the  structure  of  an 
ionizing  shock  wave.  As  discussed  in  the  previous  section,  the  Ranklne-Hugoniot 
front  is  evident  as  a sudden,  almost  discontinuous  shift  in  the  fringes.  Behind 
this  front,  changes  in  the  flow  properties  bake  place  a lowly  as  the  gas  begins  to 
relax  (this  shows  up  as  a region  in  whieh  fringe  changes  are  virtually  indiscernible). 
The  'ionization  front',  through  whloh  a rapid  increase  in  the  degree  of  ionization 
takeB  place  up  to  equilibrium  value,  is  apparent  as  a closing  together  of  the  fringes. 

It  may  be  noted  that  while  a nondimens  io.ial  fringe  shift  through  the 
Ranklne-Hugoniot  front  ia  about  0,75  for  k^  and  1.5  for  kg  for  a density  ratio  of 
about  4,  the  shifts  through  the  'ionization  front'  are  larger  (5.5  for  k^,  and  2 
for  kg)  even  thovgh  the  density  ratio  is  only  about  2.  This  may  be  aacounted  for 
by  the  negative  election  contrlbituion  to  refractive  index,  to  which  the  higher  wave- 
length is  more  sensitive. 

Analysis  of  the  interferograms  illustrated  in  Figs.  11-13  together  with 
those  of  Fig.  14,  with  the  method  described  in  Appendix  C,  yield  the  freestream 
relaxation  region  density  and  ionization  profiles  given  in  Figs*  15  and  16.  The 
values  p^,  Pg  and  p.  are  respectively  the  initial  pre-shock  gas  density,  the  frozen 
(Ranklne-Hugoniot)  density  and  the  post-relaxation  equilibrium  density.  The  limit 
of  in  the  perfect  gas  case  is  4 [i.e.jPg/Pj.  -*«»Pg/P1  -» (71+l)/(r1-l),  where 
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y is  the  specific  heat  ratio,  y * 5/3] i ?or  a monatomic  gas.  From  Fig.  15,  the 
density  ratio  P2/p^“  3.96  for  a pressure  ratio  of  P2/P1"  345  across  the  shock 

front,  while  from  Fig.  l6,  Pg/P^”  3 *9^  for  P2//pi“  2^1,  1116  real-  (i®Perfecfc)  gas 

density  ratio  p^/p^,  which  relates  equilibrium  gas  density  to  that  of  the  pre- 
shook state,  has  values  of  p /pj-  7.7**  for  Pg/Pj-  396  (argon)  and  P^/Pf  7.83  for 
Pg/P^  405  (argon  + hydrogen;  from  Fig.  15,  and  from  Fig.  16,  Pj/Pj"  6.08  for 

Pg/Pj-  258. 

A comparison  between  initial,  frozen  and  equilibrium  flow  properties  for 
pure  argon  is  made  in  Table  II. 


TABLE  lit  SHOCK  WAVS  FLOW  PROPERTIES 


Mb  - 16.53  | 

P(torr) 

p(Kg/m3) 

T(e»K) 

uOVO 

a($) 

— EM 

Initial  (1) 

5.12 

0.0111 

297 

0 

0 

0 

Frozen  (2) 

1749 

0.0438 

25,588 

3962 

0 

0 

Equilibrium  (3) 

2038 

0.0857 

13,158 

4618 

15.78 

2.039X1023 

Ma  ■ 13.59 

Initial  (l) 

5.09 

0,0110 

297 

0 

0 

0 

Frozen  (2) 

1175 

0.0433 

17,390 

3253 

0 

0 

Equilibrium  (3) 

1315 

0.0668 

11,741 

3643 

7.42 

7.47x10s2 

Figure  15  indicates  that  the  relaxation  length  for  a shook  wave  of  Mach 
number  16.53  travelling  into  pure  argon  (Pj*  5.12  torr)  is  about  2 cms.  This  is 

rather  longer  than  the  theoretical  value  predicted  by  Hoffert  and  Lien  (Ref.  1) 
who,  for  the  conditions  M ■ 16.5  and  P^*  5.0  torr,  give  a relaxation  length  of 
1.6  cms.  The  experimental  results  obtained  by  Wong  and  Bershader  (Ref.  3)  give 
a value  of  1.22  cms.  for  corresponding  conditions.  The  effects  of  radiation 
losses  in  the  post  shock-wave  equilibrium  region  are  evident  in  the  gradual  de- 
crease in  a end  p. 

The  results  for  relaxation  behind  a shock  of  Mach  number  Mg"  13.55 
moving  Into  argon  (P,«  5. 09  torr)  are  shown  in  Fig.  16.  A relaxation  length  of 
approximately  8.0  ami  is  obtained  which  compares  favourably  with  Hoffert  and 
Lien's  predicted  value  of  7.6  cm.  for  corresponding  conditions.  Wong  and  Bershader 
give  a result  of  6.1  cms  which  again  is  somewhat  lower.  The  experimental  equili- 
brium value  of  ionization  of  about  4.5$  is  considerably  lower  than  the  predicted 
value  of  7.5$.  There  is  also  a good  deal  of  scatter  in  the  expe.  lmental  results, 
notably  density.  The  usefulness  of  the  interferometric  experimental  technique 
must  be  questioned  in  flow  regions  where  density  and  degree  of  ionization  changes 
are  very  Bmall  and  the  resulting  fringe  shifts  and  refractive  index  changes 
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are  of  the  same  order  as  the  errors  associated  with  such  measurements. 


Applying  the  error  equation  given  in  Sec.  4.1  (Eq*  4,1,1)  to  the  inter- 
ferometric equations,  Eqs.  (c6)  and  (C7)  the  respective  errors  for  p and  a are 

0.5*10"**  Kg/nH  (or  2#  of  p0)  and  0.1#,  for  a refractive  index  difference  (N)  error 
-f  10“7.  From  the  theoretical  profiles  shown  in  Fig.  1$  the  average  change  in  p 
and.  a over  the  first  5 ems,  may  be  taken  as  0.0002  Kg/nH/cm  (or  0,5#  of  p„  per  cm) 
and  0,1#/ cm.  respectively.  For  an  interferogram  with  fringes  spaced  by  5.0  ran 
the  change  in  p and  a from  fringe  to  fringe  expected  is  4*10" 5 Kg/m’  (0.1#  of  p2) 
and  0.02#,  which  are  considerably  smaller  quantities  than  the  experimental  error 
values.  Consequently,  the  shorter  the  relaxation  region  (or  the  smaller  the  region 
of  small  flow  changes)  the  better  should  be  the  agreement  between  experiment  and 
theory  (with  loss  scatter  in  tile  experimental  data) » This  is  well  illustrated  by 
the  three  sets  of  results  shown  in  Figs.  15  and  l6,  the  best  agreement  being  ob- 
tained with  the  argon  + hydrogen  data,  which  will  now  be  discussed. 

Bristow  (Ref.  18)  discussed  the  use  of  hydrogen  addition  with  a view  to 
eliminating  some  inexplicable  flow  disturbances  both  in  and  behind  the  shock  wave. 
The  same  type  of  disturbances  were  observed  at  higher  Mach  numbers  in  this  series 
of  experiments.  However  since  only  a eonparatively  small  aperture  of  the  flow 
field  was  observed  in  each  experiment,  the  disturbances  were  not  so  apparent. 

Figure  17  shows  interferograma  obtained  by  Bristow  for  an  ionizing  shook  wave 
(M  ■ 17)  over  the  couple te  7 inch  aperture  of  the  shock  tube.  Some  possible  causes 
of s these  flow  disturbances  are  suggested  but  not  explained  in  Ref.  18,  The  addi- 
tion of  0.4#  Hb  (by  pressure)  to  the  argon  test  gas  was  found  by  Bristow  to  elimi- 
nate these  disturbances.  The  effect  of  hydrogen  inpurity  on  the  shook  wave  relaxa- 
tion length  is  quite  drastic  (Fig.  12).  A comparison  with  the  pure  argon  gas  case 
is  made  In  Fig.  15,  whore  it  la  seen  that  the  relaxation  length  is  reduced  from 
20  mm  to  approximately  6 mm,  for  similar  conditions.  This  may  be  explained  by  the 
larger  excitation  collision  cross-sections  for  H-H  and  H-e  as  compared  with  those 
for  A-A  and  A-e.  A full  discussion  on  the  effect  of  hydrogen  addition  is  given 
by  Liu  in  Ref.  15. 

It  is  worth  noting  that  the  equilibrium  degree  of  ionization  is  rather 
lower  (13*5#)  than  that  predicted  for  pure  argon  (16.31#)  while  density  is  higher 
(0.0975  Kg/nP  as  opposed  to  0.0873  Kg/o3  for  pure  argon).  The  major  reason  for 
this  Is  the  effect  of  the  low  hydrogen  atomic  mass  (about  1/40  of  that  for  argon), 
as  described  in  Ref.  15. 

The  theoretical  results  obtained  by  Liu  show  good  agreement  with  experi- 
ment both  for  the  pure  gas  case  and  with  hydrogen  additive. 

4,3  Wall  Effects  in  Relaxation  Region 

Perhaps  the  most  interesting  results  obtained  in  this  experimental  pro- 
gram were  obtained  in  the  relaxation  region  close  to  the  shock  tube  wall  (see 
FigB,  11-14),  These  are  believed  to  be  the  first  of  such  experimental  data  obtained 
in  this  region. 

A large  concentration  of  electrons  is  apparent  as  a sheath  growing  from 
the  shock  front,  the  degree  of  ionization  increasing  much  more  rapidly  than  in  the 
freestream.  Figures  18  and  19  show  contours  of  constant  a for  shocks  (M  ■ 16,53 
and  M ■ 13*59)  moving  into  pure  argon.  While  similar  effects  are  seen  ift  both 
instances,  the  phenomenon  is  much  more  in  evidence  with  the  stronger  shock  as  would 
be  expected,  the  relaxation  length  being  much  shorter  and  the  equilibrium  degree  of 
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ionization  being  about  twice  that  for  the  weaker  shock.  It  is  difficult  to  conclude 
from  the  interferograms  (Fig. 12)  whether  a similar  wall  effect  occurs  with  the 
addition  of  hydrogen  to  argon  test  gas.  The  relaxation  region  was  shortened  to 
the  extent  that  only  a few  fringes  transmit  data  in  the  wall  relaxation  region. 

The  fringe  spacing  should  be  reduced  considerably  from  that  shown  in  Fig.  12 
(2.24  mm  for  1.10  mm  for  \,)  in  order  to  obtain  the  required  information. 

The  reasons  for  this  premature  ionization  close  to  the  wall  are  far  from 
clear.  One  considered  possibility  was  that  a gas-surface  interaction  occurred 
between  the  argon  plasma  and  the  chromium  plated  steel  shock  tube  wall*  Two 
experiments  were  carried  out  to  try  and  eliminate  this  possiblity.  In  one,  a 
covering  (0.003")  of  ordinary  Sellopac  (cellophane)  tape  was  stuck  to  the  wall  and 
in  the  other  a thin  tungsten  foil  (0.001")  was  fixed  to  the  wall.  It  was  hoped 
that  the  tape,  being  a good  electrical  insulator,  or  the  tungsten,  having  a different 
work  function  than  chromium  or  steel,  might  change  the  electron  distribution  if 
indeed  a gas-surface  interaction  was  taking  place.  However,  no  changes  were  apparent. 
Figure  20  shows  intarferograms  of  the  shockwave  passing  over  the  tungsten  foil. 

A weak  Mach  wave  may  be  seen  emanating  from  the  leading  edge  of  the  foil  at  the 
left  side  of  the  interferogram. 

The  possibility  of  ispurities  being  present  on  the  shock  tube  wall  must 
also  be  considered.  The  wall  was  cleaned  as  thoroughly  as  possible  before  each 
experiment.  However  the  presence.of  water  vapour  molecules,  typically  at  a density 
of  1011+  cm"*2  at  a pressure  of  10’®  torr,  must  be  considered  a likely  cause,  since 
at  the  tenperaturas  behind  the  shock  front  (say  20,000°K)  any  water  vapour  molecules 
would  be  conpletely  dissociated  and  hydrogen,  as  already  proven,  has  a dramatio 
effect  on  the  ionization  rates.  However  it  is  not  clear  how  their  presence  would 
be  felt  in  the  froostream. 

It  is  of  interest  to  note  that  in  another  experimental  study,  in  the  seme 
facility,  involving  the  use  of  a finite  flat  plate  model,  a similar  phenomenon 
was  evident  (Ref.  16), 

4.4  Boundary  Layer 

Before  any  experiments  were  conducted  to  obtain  boundary  layer  data  for 
ionized  argon  gas  flows,  a check  on  the  experimental  technique  was  made  by  obtaining 
some  data  for  boundary  layers  behind  weaker,  nonionizing  shock  waves,  that  could  be 
conparcd  with  well-established  classical  theory. 

Figure  21  shows  experimental  data  and  a theoretical  density  profile  at 
10.0  cma  behind  a shock  wave  of  Mach  number  M ■ 6.48  travelling  into  argon  test 
gas  at  Initial  pressure  P ■ 3.06  mmHg.  The  experimental  results  agree  reasonably 
well  with  the  similarity  solution  obtained  by  Whitten  (Ref.  16). 

The  errors  Involved  in  these  data  concern  the  shock  tube  wall  location 
and  fringe-shift  measurements.  Determination  of  the  wall  position,  as  described  in 
Sec,  3.3.4,  is  subject  to  a + 0.05  mm  error.  Following  the  error  analysis  described 
in  Sec.  4.1,  the  resulting  error  in  density  calculations  from  fringe  shift  measure- 
ment errors  of  0.03  in  Ap  ■ + 6,6«10*4  Kg/m3,  or  Ap  ■ 1.6^  of  p , the  freestream  den- 
sity value.  Thus  the  largest  source  of  error  is  the  wall  location  error.  This  is 
representei  in  Fig.  21  by  an  error  bar. 

Figures  22  and  23  are  flow  and  no-flow  interferograms  for  A.  and  V 
respectively,  taken  for  a shock  Mach  number  of  l6.?0  and  P^«  5.63  torP.  The^flow 
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disturbances  mentioned  in  Sec.  4.2  are  very  much  in  evidence.  Care  was  taken, 
when  analyzing  the  boundary  layer,  to  avoid  regions  in  which  disturbances  were 
apparent.  Figure  24  shows  another  set  of  flow  inter ferograms .capturing  the  boun- 
dary layer  immediately  behind  a shock  wave  of  similar  strength  (M^  ■ 16.47, 

P.  = 5.22  torr) . 3 

Density  and  degree  of  ionization  data  taken  at  3.75  cma  behind  the  shock 
front  are  presented  in  Fig.  25,  together  with  the  equilibrium,  frozen  and  non- 
equilibrium  solutions  obtained  by  Liu  and  Whitten,  Details  of  their  theoretical 
solutions  are  given  in  Refs,  15  and  16,  The  experimentally  determined  boundary 
layer  thickness  (based  on  a 99$  density  ratio  value)  is  0,7  mm  while  that  predicted 
by  equilibrium  theory  is  0.9  wm.  Comparison  with  the  theoretical  models  suggests 
that  the  boundary  layer  is  either  frozen  or  in  nonequilibrium,  which  is  contrary 
to  the  as  suctions  made  by  most  authors. 

Apart  from  errors  arising  from  refraction  effects,  which  are  not  con- 
sidered in  this  study,  the  errors  in  p and  a due  to  wall  position  error  became 
quite  significant  in  dual-wavelength  interferometry  (Appendix  C),  maximizing 
towards  the  wall.  Even  for  an  error  of  only  0.0?  mm  in  wall  position,  the  result- 
ing fringe  shift  error  may  be  0.3  close  to  the  wall  giving  density  and  degree  of 
ionization  errors  of  0,01  Kg/m’  and  0,5$  respectively  or  10,7$  and  3.2$  of  their 
respective  freestresm  values. 


Even  though  refraction  errors  have  been  neglected,  the  results  of  Kuiper 
(Ref,  24),  for  fringe  shift  error  due  to  refraction  through  the  end  wall  boundary 
layer,  suggests  that  even  close  to  the  wall  (0.2  mm)  where  refraction  effects  are 
strongest,  the  theoretical  error  is  only  0,02  (for  similar  gas  conditions),  which 
is  small  compared  to  the  fringe  shift  error  due  to  wall  position  error. 

Figure  26  shows  flow  interferograas  of  the  boundary  layer  region  between 
18  and  24  cm  behind  a shock  of  Mach  number  16,49.  It  may  bo  appreciated  from  a 
qualitative  study  of  the  interf erogrwns , that  the  boundary  layer  thicker  consid- 
erably in  this  region  (to  the  left).  Measurements  indicate  that  at  x ■ 20,0  om, 
the  boundary  layer  is  approximately  2mm,  while  at  x ■ 22,5  cm  it  haB°thickened 
to  about  4 mm.  In  the  absence  of  any  thin- film  gauge  traces,  it  might  be  inferred 
from  this  rapid  thickening  of  the  boundary  layer,  that  transition  from  laminar  to 
turbulent  flow  is  taking  place. 


For  a Reynold's  number  defined  as  (Ref,  11) 


Re  ■ 


where  V2  ■ Ug-Ug,  ■ usi  a Reynold's  number  of  Re  * 2.65*10^  is  obtained. 

Experimental  and  theoretical  profiles  for  x ■ 20.0  eras  are  shovm  in 
Fig.  27  and  again  the  large  differences  between  equilibrium  theory  and  experiment 
may  be  noted. 

Additional  errors  to  those  mentioned  earlier  may  arias  in  this  case  due 
to  the  fact  that  since  the  shock  wave  is  not  in  the  field  of  view,  the  exact  de- 
termination of  x is  not  possible.  However  from  recorded  time  measurements,  i.e,, 
the  times  betweeh  the  shock  wave  passing  through  the  centre  of  the  test  section  and 
and  the  time  at  which  the  Pockels  cell  is  triggered,  the  position  of  the  Bhock  wave 
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with  respect  to  the  field  of  view  may  be  determined  to  within  + 0.5  cms  or  better. 
The  errors  incurred  in  p and  a are  virtually  negligible. 

Figure  28  shows  results  taken  12  cms.  behind  a shock  wave  of  M “ 13.56, 
where  the  experimental  density  boundary  layer  thickness  is  about  1.2  mm  while 
a boundary  layer  thickness  based  upon  a degree  of  ionization  ratio  a/<>  “0.99 

appears  to  be  significantly  greater  at  approximately  1.6  mm.  This  sugfests  the 
presenc  - of  a thermal  boundary  layer  which  is  rather  thicker  than  the  viscous 
boundary  layer. 

5.  CONCLUSIONS 


The  use  of  duel -wavelength  interferometry  has  proven  successful  in  the 
study  of  ionizing  shock-wave  structure  and  the  shock-tube  boundary  layer  in  argon. 

Results  taken  of  the  shock  front  thickness  indicate  that  the  inter- 
ferometer-light source  alignment  procedure,  described  in  Appendix  B,  should  be 
extended  to  improve  alignment  about  the  y-axis  of  the  shock  tube. 

Although  the  effects  of  refraction  through  the  shock  front  are  considered 
to  be  slight  for  a perfectly  aligned  interferometer,  a full  analyses  should  be 
carried  out  to  determine  tho  magnitude  of  refraction  errors  for  a misaligned 
interferometer. 

The  measurements  of  density  and  ionization  in  the  relaxation  region  show 
good  agreement  with  theory  although  the  scatter  in  experimental  results  for  weaker 
shocks  (M  m 13.5)  suggest  that  further  work  should  be  directed  at  stronger  shock 
studies  o?  that  the  laser  ruby  element  (69^3*)  be  replaced  by  & neodymium-glass 
element  (10600  X)  for  greater  sensitivity  to  electron  density  for  weaker  shocks. 

The  contribution  of  excited  state  refract! vi ties  for  neutral  argon  and 
the  effect  on  p and  a calculations  should  be  investigated  for  strong  shocks  (Ref.l6\ 

The  phenomenon  of  a faster  ionization  rate  taking  place  close  to  the 
wall  was  an  unexpected  result  of  this  experimental  investigation.  A more  detailed 
study  of  this  should  be  carried  out  to  establish  whether  the  phenomenon  is  peculiar 
to  the  particularly  larger-area  shock  tube  used  in  this  experimental  program  or 
whether  it  applies  to  shock  tubes  of  smaller  size,  where  this  phenomenon  has  appa- 
rently not  teen  noticed  (e.g. , Ref.  3). 

The  results  obtained  for  the  induced  boundary  layer,  when  compared  with 
the  equilibrium,  frozen  and  nonequilibrium  theoretical  profiles  of  Liu  and  Whitten, 
suggest  that  the  boundary  layer  is  in  a frozen  or  nonequilibrium  state.  The  latter 
profiles  are  very  similar  and  it  is  virtually  impossible  to  distinguish  between 
them  experimentally. 

Although  the  experimental  results  presented  take  no  account  of  refraction 
effects  in  the  boundary  layer,  the  theoretical  error  results  presented  in  Ref.  2l 
indicate  that  errors  due  to  refraction  are  small  compared  to  the  wall  location  error. 

The  freestream  conditions  used  in  the  calculations  were  those  predicted 
by  equilibrium  strong-shock  theory  (Ref.  16).  Any  losses  due  to  radiation  in  the 
freestream  were  neglected.  Results  presented  in  Ref.  25  for  a shock  of  Mach  number 
l6.3  and  P.0  5 oorr,  indicate  that  the  degree  of  ionization  drops  from  an  equilibrium 
value  of  about  15$  immediately  behind  the  shock  wave  to  10$  at  5 cm  behind  the  shock, 
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after  which  it  remains  fairly  constant.  Corresponding  changes  in  density  are 
relatively  small  (typically  % of  the  equilibrium  value,  p,,).  Reference  22  gives 
some  data  for  radiation  losses  obtained  in  the  facility  used  for  this  study, 

pertaining  to  stronger  shocks. 

Although  tentative  indications  of  transition  in  the  boundary  layer  were 
present  for  the  stronger  shock  case,  in  the  form  of  a rapid  growth,  no  definite 
proof  was  obtained.  A further  study  of  transition  behind  ionizing  shocks  should 
include  the  use  of  a thin- film  heat- transfer  gauge. 

The  shock-transition  density  and  degree  of  ionization  profiles  as  well 
as  those  through  the  boundary  layer,  together  with  similar  contours  of  the  rapid 
ionization  region  near  the  wall  should  provide  useful  data  for  additional  theoreti- 
cal research. 
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APPENDIX  A:  THERMODYNAMIC  AND  TRANSPORT  PROPERTIES  OF  PARTIALLY  IONIZED  ARGON 


i 


One  of  the  major  difficulties  in  solving  the  boundary  layer  equations 
for  an  equilibrium  partially  ionized  gas  flow  in  the  shock  tube,  is  the  determi- 
nation of  the  thermodynamio  and  transport  properties  (especially  the  latter)  over 
the  wide  range  of  tender atures  encountered  through  the  boundary  layer,  from  the 
aold  (300°K)  wall  to  the  hot  (typically  1? ,000°K)  free-stream. 

Some  experimental  data  have  bean  obtained  (e.g.,  Ref.  Al)  in  the  shock 
tube,  utilising  the  hot  slug  of  stagnant  gas  behind  the  reflected  shook  wave. 
However  the  majority  of  information  provided  in  the  literature  is  of  n theoretical 
nature,  sufficient  to  predict  the  thermal  conductivity  and  viscous  properties  with 
a reasonable  degree  of  accuracy. 

Previous  assumptions  such  as  a constant  Prandtl  number  and  density- 
viscosity  product  are  now  known  to  be  inaccurate, 

The  degree  of  ionization  in  equilibrium  flow,  for  given  pressure  and 
temperature,  is  given  by  the  Saha  equation  as 


p n 

where  p is  in  N/m  and  the  characteristic  temperature  of  ionization  is  182,900  K. 

Following  Ref.  A2,  enthalpy  may  be  written  as  (neglecting  excited  state 
contributions) 


h - R [5/2  (l  + a)  T + a Sj  ] 


(A2) 


and  the  equilibrium  specific  heat  at  constant  pressure, 

dh  i 


as 


Op  - * [s/s  <1  + o)  + S/a  * (§)  4 4 (*)  ] (A3) 


From  (Al), 


ajj  l - 1.25  p t*7/2  e 1 (2.5  + ^ ) 


(A 4) 


and  substitution  of  (A4)  into  (A3)  gives 


Al 


(A5) 


°P  " R { 5//2  ^ + + (5/2  T + 9j) 

[l.25  O3?  T‘7/2e  1 (2.5  + J } 

The  collision  cross-section  date  required  for  the  calculation  of  viscosity 
and  thermal  conductivity  may  be  sunnarited  as  follows t 

i)  Atom-atom  collision  cross-section  is 

SAA  " 17  ( ~I?)  °,26x  10“18  ft2  (A6) 

from  the  curve  fitted  (Ref.  A2)  to  the  data  provided  by  Ref.  A3 

11)  The  Coulomb  scattering  crop e-section,  is  (from  Ref.  Ah) 


where  m is  the  reduced  mats  (in  this  case  mA/2),  v is  the  random  velocity  given 
by  r a r 

'AHffT 

and  n is  electron  number  density, 

i*8‘*  q P 

ne  " TTa  » 

q is  electron  charge  (1.6021  lO*1^  Coulomb),  € is  permittivity  of  vacuum 
(8.851*  IQ-12  farad/m)  and  k,  the  Boltamann  constant  (l.380t  lO'2^  joule/°K). 
Substitution  of  the  constants  in  Eq.(A7)  yields 

S,  ■ h!!L^  *»  [e  W (AB) 

where  the  electron  number  density  may  be  calculated  from  Eq.  (Al). 

ill)  The  atom-ion  hard  Bphere  collision  cross-section  is  found  to 
be  virtually  independent  of  temperature  and  a value  of 

■ 1.3 5x10' 18  m2  (A9) 


A2 


y 


s 


is  used  here. 

Viscosity,  |i  may  then  b*  written  as 


l + 

a 

T3r 

V 

«AA 

1 + 

0 

TS 

O Su, 


SI  “o 

32  *AA 


following  Ref.  A2*  Any  contributions  made  to  viscosity  by  the  electrons  are 
neglected* 

In  Eq.(AlO)  the  mean  thermal  atonic  speed  is  given  by 

n . f LB  f/2 
A L » m,  J 

The  total  thermal  conductivity  K-—,  ,1s  given  by  the  mixture  rule 


where  K , the  thermal  conductivity  of  the  fully  ionized  gae  is 
s 

K.  ■ 3.55«10'17  T n*/2  cal/em  sec°K 

8 6 

The  electron-atom  average  momentum  transfer  cross-section,  Q . may  be 
obtained  by  curve  fitting  data  given  by  Ref.  a6. 
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APPENDIX  Bj  ALIONMENI  OF  INTfflFEROMBm  AND  LIGHT  SOURCE 


With  the  use  of  interferometry  for  the  atudy  of  boundary  layers,  the 
importance  of  precise  alignment  of  the  instrument  and  its  light  source,  with 
respect  to  the  surface  of  interest,  is  critical.  Any  misalignment  of  the  light 
beam  with  respect  to  the  surface  results  in  the  formation  of  spurious  fringe 
shifts  at  the  wall,  which  interfere  with  accurate  evaluation  of  the  inter  ferogram,, 
in  the  boundary- layer  region. 

The  Mach-Zehnder  interferometer  used  in  this  experimental  study  is  well 
described  In  Ref.  Bl.  It  is  a relatively  large  piece  of  equipment  (as  can  be  seen 
from  Fig.  3) » having  a 9- inch  diameter  field  of  view.  Consequently  some  care  is 
required  in  obtaining  correct  adjustment  of  It  with  respect  to  the  shook  tube  test 
section. 


The  procedure  adopted  is  described  here  step  by  atap  (with  reference  bo 
Figs . 5 i»  9) • A gas  laser  is  used  to  provide  a low  intensity  monochromatic  light 
source,  input  to  the  interferometer  is  via  a plane  mirror 

1)  Both  the  laser  end  interferometer  are  crudely  adjusted  so  that  the 
laser  beam  passes  approximately  through  the  centre  of  the  elements  of  the  inter- 
ferometer end  through  the  oentre  of  the  area  of  Interest  in  the  test  section, 
i.e.,  Just  above  the  shook  tube  floor  end  midway  along  the  visible  portion  of  it, 

ii)  Cross-hairs  are  fixed  on  the  first  parabolic  mirror  P.ao  that  the 
laser  beam  la  centered  on  them.  A eat  of  croaa-hairs  is  then  positioned  on  the 
second  parabolic  mirror  ?9,  centered  at  the  corresponding  point  on  the  mirror  as 
the  first  set. 

ill)  The  focal  point  cf  P,  is  looated  using  a portable  white  light  source 
clipped  to  the  front  of  the  iris,  17  (which  is  attached  to  the  interferometer), 
the  reflections  of  the  point  souroe^from  the  test  section  end  compensating  chamber 
windows  being  brought  to  a focua  on  the  iris.  The  positions  of  the  reflections 
with  respect  to  the  centre  of  the  iris  are  noted  and  the  following  steps  are  taken* 

lv)  The  gas  laser  is  now  adjusted  so  that  the  beam  passes  through  the 
oentre  of  the  Iris  onto  the  cross-hairs  of  F^. 

v)  P,  is  adjusted  to  direct  the  beam  onto  the  cross-heirs  of  P_  (which 
results  in  the  Changing  of  position  of  the  focal  point  of  P.  since  the  incident 
beam  to  P^  ia  off-axis),  1 

Steps  (ill)  to  (v)  are  repeated  until  the  reflections  on  I.  are  centered 
as  closely  as  possible. 

At  this  point  one  should  ensure  that  the  beam  is  still  centered  through 
the  required  part  of  the  test  section.  It  may  be  necessary  to  move  the  inter- 
ferometer and  light  source  laterally  and  vertically  to  re-centre  the  beam  through 
the  test  section. 

Assuming  that  the  Interferometer  is  in  adjustment  for  producing  white 
light  fringes  (the  procedure  for  which  is  described  in  Ref.  Bl),  a strong  white 
light  point  source  ia  placed  at  the  position  of  I.  and  the  fringes  focus  Bed  on  a 
screen  placed  at  the  image  focus  plane  of  Pg,  as  A described  in  Sec, 


Bl 


The  white  light  source  is  then  replaced  by  the  gas  laser  and  & lens 
positioned  to  focus  the  beam  to  a point  at  X^. 

Spurious  fringes  at  the  wall  indicate  that  the  interferometer  is 
mis-aligned  in  pitch.  This  may  be  corrected  by  adjustment  of  the  jacks  or 
the  tire  pressures  of  the  pneumatic  tires  on  which  the  interferometer  stands. 

Any  adjustment  of  the  interferometer  in  this  manner  requires  that  the  light 
source  be  re-adjusted  so  that  the  laser  beam  passes  through  the  iris  centre  onto 
the  first  set  of  cross-hairs. 

The  camera  is  now  positioned  so  that  the  beam,  which  exits  the  inter- 
ferometer via  the  plane  mirror  Mj,»  passes  through  the  centre  of  its  shutter 
and  the  components  of  the  cameraHare  positioned  accordingly,  so  that  the  beam 
is  directed  to  the  centre  of  the  film  plate  holders.  Using  a white  light  source 
the  camera  is  focussed  in  the  desired  plane  (see  Sec.  3,3,4)  and  the  fringes 
focussed  in  that  plane  by  remote  control  of  the  motor  driver  splitters  (3  and 
Sg)  end  the  plane  mirror  M^.  1 

The  pulsed  laser  is  aligned  by  means  of  cross-hairs  attached  to  both 
ends  of  it,  they  being  centered  on  the  optical  axiB  of  the  ruby  lasing  element, 

Polaroid  snapshots  are  taken  (Pig.  Bl)  to  ensure  that  a satisfactory 
alignment  has  been  achieved.  The  3471, 5a  component  fringes  are  focus b 3d  in  the 
camera  by  a trial  and  error  method  of  moving  the  film  holder  and  taking  several 
Polaroid  snapshots  until  the  sharpest  focus  of  the  'Invisible'  fringes  is  obtained. 
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APPENDIX  Ct  ANALYSIS  C?  IlttERFEROORAMB 


The  task  of  reducing  an  interferogram,  which  basically  consists  of  line: 
of  constant  order  of  interference,  to  useful  information  about  particle  density 
(atom,  ion  or  electron),  is  a somewhat  tedious  one,  but  accuracy  e.t  this  stage  i 

essential  to  tha  overall  success  of  an  ejtperimont. 

The  degree  of  difficulty  and  the  amount  of  error  which  is  acceptable, 
varies  with  the  particular  region  of  flow  being  explored  and  the  emount  of  time 
available  for  the  data  reduction  prooess.  In  thia  experimental  program  two  basic 
regions  were  investigated,  i)  the  ehock  wave  structure,  both  in  the  free-stream 
and  in  the  neighbourhood  of  the  wall  and  li)  the  equilibrium  boundary  layer  region 
behind  the  shock  wave.  For  the  former  case,  which  is  essentially  a 2- dimensional 
problem,  when  wall  effects  are  eonaidered,  density  gradients  are  experienced  both 
parallelto  the  flow  direction  and  perpendicular  to  it.  Whereat  in  the  equilibrium 
boundary  layer  region  where  one  assumes  constant  free-streem  flow  conditions,  we 
have  at  any  given  station,  x,  a one- dimensional  problem  with  density  changes  per- 
pendicular to  the  wall  along  y.  Consequently,  it  was  found  to  be  convenient  to 
deploy  two  Methods  for  analysing  the  interferograms,  which  will  be  described  here. 

For  parallel  fringe  Interferometry,  the  equation  relating  refractive 
index  n at  a point  (x,y)  to  fringe  order  N is  (Ref,  Cl) 

n(x>y)  £ + y [ (f  + const - n^J  ■ N (x,y)  (Cl) 


where  L is  the  test  section  width,  X is  wavelength  of  light,  d la  fringe  spacing 
(taken  from  the  no-flow  interferogram)  and  tho  suffix  'O'  denotes  reference  region 
values . 


Equation  (Cl)  is  valid  for  any  general  point  (x,y)  in  the  flow  field. 
Considering  now  e reference  point  (xQ,yo)  then  we  may  write 

n(xQ,y0)  L + [ (a2  + oonot')  - v]"  N(*o‘yo*X  (C8) 

Subtracting  (C2)  from(Cl)  gives 

[n(x,y)  - n(xQ,y0)j  j ■ N(x,y)  - N(x0,yQ)  - (C3) 

At  this  point  it  is  pertinent  to  dlscuse  the  initial  orientation  with 
which  the  fringes  are  pre-set.  The  two  basic  choioes  are  to  have  the  fringes  set 
parallel  to  the  direction  of  the  density  change  or  perpendicular  to  it.  The  former 
gives  the  moet  information  since  one  has  a continuous  trace  of  refractive  index  ir 
the  direction  of  the  changes  while  the  latter  gives  only  discrete  points  but  is 
much  simpler  to  evaluate.  In  the  shock  wave  region  there  is  no  well  defined  point 
or  region  at  which  the  flow  conditions  are  known,  as  required  by  Eq.(C3).  One 
could  assume  the  froaen  flow  values  predicted  by  theory  for  the  region  ismediately 
behind  the  shock  front.  However  the  pre-ihock  gas  state  is  well  known  from  labo- 
ratory measurements,  so  this  may  be  used  as  a good  reference  region.  In  order  bo 


Cl 


use  the  undisturbed  gas  as  a reference,  the  fringes  should  be  orientated  so  as 
to  pass  through  the  shock  front,  making  horizontal  fringes  most  appropriate.  How- 
ever, in  the  boundary  layer  region  at  the  wall  it  is  desirable  to  have  the  fringe 
orientated  vertically.  Consequently  a coinpromise  between  these  two  requirements 
is  best  obtained  by  using  diagonally  orientated  fringes. 

Looking  now  at  Eq.(C3)  in  context  with  Fig.  Cl  we  see  that  the  refractive 
index  n at  point  (x,y)  in  the  flow  field  may  be  obtained  by  the  following  pro- 
cedures 

i)  Select  a fringe  in  the  no- flow  field  to  use  as  a reference.  Arbi- 
trarily set  the  fringe  order  (e.g.,  N ■ 1)  since  it  is  AN  that  we  are  interested 
in. 


ii)  For  y ■ y , measure  x-xQ  and  calculate  AN  (in  this  case,  l),  being 
aareful  with  the  sign  of  AN. 

iii)  Measure  d,  the  no-flow  fringe  spacing  (in  x co-ord.  direction). 

iv)  Substitute  the  above  values  in  Eq.  C3  to  obtain  n(x,y)  - n(x0,yQ), 

(y  - yD). 

The  relationship  between  refractive  index  and  density  is  (Ref.  CS) 

n - 1 - Ka  p0  «*> 

where  K is  the  Oladstone-Dale  constant  for  species  's'  and  Pg  is  the  density  of 
speoies  's'«  When  expanded  for  the  case  of  a singly  ionised  plasma,  neglecting 
contributions  to  refractive  index  from  excited  states,  Eq.  C4  becomes  (in  e.g.s. 
units) 

(nJ-l)-(nl-l)  ■ Ka  (Pj"Pi)-(Pjaj-P101)(KA-Ki  * 0*67«109  X2)  (C5) 

where  'i'  denotes  reference  value,  'j'  the  unknown  value  and  and  are  the 
Glads tone-Dale  constants  for  atoms  and  Ions  respectively. 

From  Eq.C5  we  see  the  requirement  for  two  similar  equations  (one  each 
for  \j_  \)  to  be  solved  simultaneously  for  pj  and  ay 

If  we  let 

N.  • [(n.  -1)  - (n.  - 1)] 
l J1 

and  Ng  - [(n^-l)  ’ (n^-  D3 

then  algebraic  manipulation  of  the  two  equations  (C5)  yields 


C2 


and 


a 


J 


(C7) 


where 


0.33  Ka 
0.33  Ka 


cl% 


+ 0.67*109 

+ 0.67**109 


* 


And  Kj  ■ 0.67  Ka  (Refs,  C2  And  C3). 

For  the  light-source  wavelengths  used  in  these  experiments 
(X^  ■ 6943  a,  ^ " 3*71.5  a)  the  constants  take  the  following  values 


K4  ■ 0.15726  cm3/gm  K.  ■ 0.16273  cm3/gm 

1 a e 

c^  ■ 3.2616  onr/gm  cg  ■ 0.8611  cm-ygm 

K « 0.3986  cm6/gm2 


Hence  with  the  values  of  N,  and  N-  taken  from  the  two  interferograms 
using  Eq.C3  one  can  obtain  p and  a from  Eqs.  C6  and  C7. 

For  the  equilibrium  flow  boundary  layer  region  the  following  procedure 
was  adopted. 

Using  fringes  orientated  perpendicular  to  the  wall,  the  change  in 
density  through  the  boundary  layer  is  characterized  by  a shifting  or  bending 
of  the  fringes  to  the  left  or  right  (depending  on  how  the  interferometer  is 
adjusted).  Teklng  the  free  stream  flow  as  our  reference,  the  equilibrium 
conditions  being  predicted  by  strong- shock  real  gas  theory,  a no-flow  interfero- 
gram  is  used  to  obtain  the  undisturbed  gas  fringe  path  to  the  wall  and  then 
using  this  in  conjunction  with  the  flow  interferogram,  the  fringe  shift  may  be 
determined  at  any  point  on  a fringe. 

The  relationship  between  the  non-dimensional  fringe  shift  S.,  (non- 
dimansionalized  by  the  no-flow  fringe  spacing)  and  ahange  in  refractive  index 
is  (Ref.  Cl) 


Su  ■ 5 C(>yi)  - (v1}  ] (C8) 


Referring  to  Fig.  C2,  the  nondlmensional  shift  S.u  is  determined  in 
the  following  manner  j 


C3 


i)  For  the  point  of  interest  (x,y)  the  average  x co-ordinate  difference  in 
the  uniform  flow  region  between  the  flow  fringe  on  which  the  point  lies  and  the 
nearest  no-flow  fringe,  is  determined  (Ax). 

ii)  At  the  same  y co-ordinate  (y  y ) the  x co-ordinate  difference  between 
the  flow  and  no-flow  fringes  in  the  boundary  layer  is  measured  (x-xQ) . 

iii)  The  non-dimensional  shift  is  then  calculated  as 


x-x  + Ax 


for  the  point  (x,y),  where  d is  the  average  no-flow  fringe  spacing. 

An  equation  similar  in  form  to  Eq.  C5  may  be  obtained  as 


8ij  ■ * iV'r'i'-'VrVi1  <v*i  ♦ °-S7,1°9 


(in  c.g.s.  unite) 

with  corresponding  equations  for  p and  a being 

. . . \C2S1  A VlS2 

Pi  " Pi  ’ 


r \\h  / 

" [ piai  ■ “Tar"  4 “hre — J / p, 


(C9) 


(CIO) 

(oil) 


where  L ia  the  test  section  width  taken  as  L * 10.16  cm,  the  other  constants  taking 
on  the  same  values  as  for  Eqs.  C6  and  C7. 

Thus  having  obtained  values  of  S,  and  S2  at  a point  (x,y),  the  densities 
(p  and  a)  at  that  point  may  be  calculated  from  Eqfi.  CIO  and  C1L. 

Some  points  should  now  be  made  concerning  the  practical  application  of 
the  described  methods  of  analysis. 

The  most  important  is  that  the  values  of  N.  and  Ng  or  S,  and  „g  may  only 
be  determined  from  the  lnterferograms  at  points  (x,y)  which  actually  lie  on  a 
flow  fringe  sinoe  it  is  convenient  to  use  fringe  oentres  as  references  of  constant 
intensity.  The  fringe  pattern  for  each  wavelength  is  of  course  different  and 
the  probability  of  a point  (x,y)  chosen  for  a fringe  at  one  wavelength,  lying  on 
a fringe  in  the  other  wavelength,  is  very  small.  Consequently  an  interpolation 
routine  ia  required  to  determine  the  values  of  and  Ng  or  and  Sg  at  points 
between  fringes.  £ 

A Ruscom  digitizer  (Model  21),  measuring  in  hundredths  of  an  inch,  was 
used  to  obtain  (x,y)  co-ordinate  data  for  each  fringe  of  an  interferogram.  Thin 
data  (punched  onto  computer  cards)  was  then  used  with  a computer  program  written 
(for  the  UTIAS  IBM  1130  computer)  to  (see  Fig.  C3) » 
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i)  Linearly  interpolate  the  digitized  x,y  co-ordinates  in  ateps  of  y 
(Ay  * 0.05  nun  in  lab.  co-ord). 

ii)  Confute  and  Ng  or  and  Sg  at  interpolated  x,y  co-ordinates. 

iii)  Linearly  interpolate  between  two  nearest  data  points  lying  either 
side  of  required  x station,  at  each  step  value  of  y for  each  wavelength. 

iv)  Calculate  p and  a from  interpolated  values  of  N.  and  Np  or  31  and 

Sg  at  each  point  (x,y).  1 2 

The  importance  of  correctly  determining  the  true  wall  position  (as 
mentioned  in  Sec.  3*3»*0  to  the  accuracy  of  p and  a,  may  he  appreciated  here. 

Any  error  in  the  wall  location  for  one  wavelength  interferogram  results  in  the 
values  of  N or  S being  assigned  the  wrong  (x,y)  values  for  that  interferogram. 
Estimates  of  the  errors  involved  in  calculated  values  of  p and  ct  due  to  wall  loca- 
tion errors  are  given  in  Sec.  U.3. 

One  basic  difference  between  the  methods  employed  is  the  use  of  the  no- 
flow interferograms , taken  for  each  experiment  prior  to  the  run.  The  purpose  of 
these  interferograms  is  to  eliminate  any  fringe  distortions  due  to  optical  imper- 
fections of  the  interferometric  system.  For  the  method  used  to  evaluate  tlie 
equilibrium  boundary  layer  region,  the  no-flow  interferogram  waE  employed  as 
previously  described.  However,  for  the  shock-wave  region  it  was  found  impractical 
to  use  the  no-flow  interferograms  owing  to  the  gross  changes  of  the  fringes  in 
the  field  of  view.  But  since  the  reference  region  taken  in  this  method  was  a 
no-flow  region,  account  was  taken  of  any  aberration  (lens)  effects  but  not  of  any 
local  fringe  distortions  due  to  burns  on  the  optical  elements  or  diffraction  caused 
by  "oreign  particles.  This  is  felt  to  be  a minor  source  of  error  and  not  signi- 
ficant to  the  experimental!  results  in  a region  where  relatively  large  changes  in 
flow  properties  take  place. 
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FIG.  21  BOUNDARY  LAYER  DENSITY 
PROFILE  K).0  CMS.  BEHIND  SHOCK 
Mg  * 6.48 


FLOW  CONDITIONS 
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FIG.  25  BOUNDARY  LAYER  DENSfTY  AND  DEGREE  OF 
IONIZATION  PROFILES  3.75  CMS.  BEHIND  SHOCK  t 
Mt- 16.50 


FLOW  CONDITIONS 
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